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Multirate filter as well as display system and mobile telephone comprising said multirate 
filter. 



The invention relates to a multirate filter as well as a display system and a 
mobile phone comprising a multirate filter. 

Digital filters find widespread use in audio and video processing, for example 
in mobile phones, set top boxes, digital television sets, and other consumer or professional 
5 products. Symmetric filters form an important class, because of their linear phase property 
and the possibility to exploit this symmetry to simplify the architecture of the filter and 
therewith reduce costs. 

In particular, multirate filters are used in applications where the output signal 
and the input signal of the filter should have mutually different sample rates. Such filters are 
10 applied for example in image processing to effect a scaling of a digitally encoded image. One 
of the most important concepts in multirate filtering is the polyphase decomposition and the 
closely related polyphase structure. This concept allows for very efficient implementations 
both in hardware and software of interpolating and decimating filters. 

A straightforward application of the polyphase decomposition however 
15 usually introduces asymmetric polyphase components, which substantially reduces the 
efficiency of the implementation 

It is a purpose of the invention to provide an architecture for a multirate filter 
wherein the symmetry of its components is recovered. For that purpose the multirate filter 
according to the invention has a construction as defined in claim 1 . It has been recognized by 
20 the inventors that a symmetric multirate filter may be constructed from an input unit, a filter 
unit and an output unit, wherein the filter unit has symmetric modules derived from the 
polyphase components of the multirate filter provided that the modules are provided in pairs 
having transfer functions Ho(z) and Hi(z), which are derived from a basic transfer function 
Hb(z) as follows. The first one Ho(z) of the transfer functions is based on the sum of the basic 
25 transfer function H B (z) and its mirrored version: 
H 0 (z) = c 0 (H B (z) +M ay ,H B (z)) , and 

the second one Hi(z) is based on the difference of the basic transfer function H B (z) and its 
mirrored version: 
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H x (z) = c x (H B (z) -M fl / a (z)) . 

The mirror operation Mo, v on the basic transfer function H B (z) is defined as: 
M ajr tfl b 00) = ccb-**H b * (z" 1 ) , and wherein 

= £V[ m l z ~ m • being the z-traraform of A d [m] 

Therein the value 2vj/ is an integer number, selected from Z, and 
the value a is an element from the set C of complex numbers. 

More in particular it has been discovered by the inventors that the asymmetric components of 
the polyphase filter can be redesigned into one of the embodiments described in claims 2 and 
5, or a combination thereof. 

In the embodiment of claim 2 the basic transfer function H B (z) is a polyphase 
component of the multirate filter. The embodiment of claim 2 has two sub-embodiments as 
defined in claim 3 and claim 4. In the subembodiment of claim 3 the input unit comprises a 
combination unit as defined in claim 1, and claim 4 the output unit comprises such a 
combination unit. In the embodiment of claim 5 the basic transfer function H B (z) is a linear 
function of two polyphase components is of the multirate filter. 

The invention further relates to a method of designing a multirate filter 
according to the invention. Such a method could either be part of a design tool, but could 
alternatively be part of a device having an adaptive filter. In that case the method enables that 
device to replace the polyphase components which it has calculated by and equivalent set of 
symmetrical modules. 

These and other aspects of the invention are described in more detail with 
reference to the drawings. Therein: 

Figure 1 illustrates a display system comprising a multirate filter according to 

the invention, 

Figure 2 illustrates a mobile phone comprising a multirate filter according to 

the invention, 

Figure 3 A illustrates a filter in the form of a tapped delay line, 
Figure 3B illustrates a simplified version of the filter of Figure 3 A, 
Figure 3C illustrates a filter in the form of an adding delay line, 
Figure 3D illustrates a simplified version of the filter of Figure 3C, 
Figure 4A shows a decimating filter, 
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Figure 4B shows implementation of the filter of Figure 4A in the form of a 

polyphase structure, 

Figure 5 A shows an interpolating filter, 

Figure 5B shows implementation of the filter of Figure 5 A in the form of a 

polyphase structure, 

Figure 6A schematically shows a multirate filter according to the invention, 
Figure 6B shows a first embodiment of the multirate filter according to the 



invention, 
invention, 
invention, 



Figure 6C shows a second embodiment of the multirate filter according to the 

Figure 6D shows a third embodiment of the multirate filter according to the 

Figure 7 A shows a filter having a rational decimation factor, 

Figure 7B shows an implementation of this filter according to the invention. 



FIG. 1 illustrates an example display system a having an image processor b 
that includes a multirate filter e. The filter e is implemented as a polyphase filter and allows 
both interpolation and decimation. Typically, the display system a includes filters for 
effecting both horizontal and vertical scaling. The controller c determines the appropriate 
scaling in each dimension, sets an appropriate mode of the filter e (interpolate or decimate, 
discrete phase or continuous phase, as required), and determines the appropriate coefficients 
that are provided by the memory d, depending upon the scaling and the mode. 

Figure 2 shows another application of a multirate filter according to the 
invention, in this case as filter for eliminating noise and other interference from a received 
signal in a mobile telephone h. A signal carrying, for example, encoded speech is received by 
antenna g, downconverted to a baseband and amplified, if necessary, by receiver h then 
routed through an analog-to-digital converter (ADC) i for conversion to discrete samples. 
The discrete samples are routed through the multirate filter j which filters out frequencies 
below a selected threshold. The threshold is selected to distinguish between information 
components of the transmitted signal and noise and other interference components. The 
signals are then demodulated by a demodulator k and decoded by a signal decoder 1. Output 
signals from the signal decoder I are routed to a speech decoder m for conversion to digitized 
voice signals. The digitized voice signals are converted to analog signals by a digital to 
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analog converter (DAC) n for ultimate output through a speaker o of the mobile telephone. 
Other components, such as error detection and correction components, may be provided as 
well within the receive portion of the mobile telephone. 

Likewise the multirate filter according to the invention is applicable in the 
sending part of a mobile telephone. The architecture of the sending part (not shown) of a 
mobile telephone is globally the inverse of the receiving part. Hence, such a mobile telephone 
comprises an analog to digital converter for converting an analog speech signal into a digital 
speech signal. A speech encoder subsequently compresses the digital speech signal and 
provides the compressed signal. A signal encoder performs a channel encoding operation at 
the compressed signal and provides a channel signal. A modulator modulates the channel 
signal. The channel signal is filtered by a multirate filter according to the invention as 
described in more detail below. The filtered signal is then converted into an anlog signal by a 
digital to analog converter. A transmitter transmits the analog signal. 

In the present description a filter h[n] having z-transform H(z) is defined as 
(a,\|0-symmetric if it fulfills the following relationship: 
(5) H (z) = az^H* (z" 1 ) , which corresponds to 

h[n] = ah* [2y/ - n] in the time domain 

A special case is (<x=l,i|/), h[n] being a real-valued function. For these filters 
h[n] the following relation holds: 
h[2y/ - i] = h[i] wherein i is an integer, 

Examples are h[]= <a,b,c,b,a> having symmetry (ot=l,i|/=2), and 
h[]= <a,b,c,c,b,a> having symmetry (<x=l,\|/=5/2) 

Another special case is (a— l,v|/). For these filters either 
h[2y/ - i] = -h[i] wherein i is an integer, 

These filters are also referred to as anti-symmetric filters, not to be confused 
with asymmetric filters. 
Examples are: 

h[] = <a,b,0,-b,-a> having symmetry (a— l,\j/=2), and 
h[] = <a,b,0,0,-b,-a> having symmetry (a— 1,\|/=5/2) 

The sign of symmetry is however not limited to real values, as is illustrated by 
the following examples: 

The filter h[] = <a+jb,b+ja> having symmetry (a=j,\|/=l/2), and 
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h[]= <a+jb,-b-ja> having symmetry (a— j,vp=l/2). 

Filters having (a,\|/)-symmetry as defined above can be implemented 
efficiently. It is known to exploit their symmetry by reducing the number of components. 
This illustrated in Figures 3 A-3D. Figure 3A shows a symmetric filter implemented in the 
form of a tapped delay line. In the example shown, the delay line comprises 4 delay elements 
z x 9 a first multiplier having multiplication factor a is coupled to the input of the delay line, a 
second, a third and a fourth multiplier for multiplying with multiplication factors b,c,d are 
coupled to subsequent connection nodes between the delay elements z" 1 , and a fifth multiplier 
having multiplication factor e is coupled to the output of the delay line. The sum of the output 
signals of the five multipliers is the output signal of the filter. The filter of Figure 3 A can be 
reduced to the filter of Figure 3B having only three multipliers if the multiplication 
coefficients a and e are equal on the one hand and the multiplication coefficients b and d are 
equal on the other hand. 

Analogously Figure 3D shows an efficient implementation of the adding delay 
line of Figure 3C. 

Another known technique in filter design is the decomposition into polyphase 
components. This is in particular relevant for decimating or interpolating filters, or filters 
using a combination of decimation and interpolation. 

Any filter H(z) can be decomposed in R polyphase components H R:r (z), which are related to 
H(z) as: 

R-\ 

H(z) = ]T z~ r H Rj . (z*) , wherein the polyphase components relate to the transfer function h() 
as: 

n 

Such a polyphase decomposition of the filter has the advantage that the 
polyphase filter components can operate at a lower clock speed than a filter which is not 
decomposed. 

This is illustrated in Figure 4A and 4B. Figure 4A shows a decimating filter 
H(z) with decimation factor D. The filter has to operate at the same clock speed as the signal 
X(z) processed by the filter. Figure 4B shows the decomposition of this filter into polyphase 

components H D :o(z), H D: i(z), , H D: d-i(z). The polyphase components are coupled via a 

decimator to a delay line having delay elements with delay z" 1 . The decimators decimate the 
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sampling frequency of the signals obtained from the delay line with a factor D, so that the 
polyphase components may operate now at a corresponding lower frequency. 

Figure 5 A and 5B illustrate a second example. In this case the output signal 
Y(z) is obtained from the input signal X(z) by an interpolating filter H(z) with interpolation 
factor I. In Figure 5 A, the interpolated signal received by the filter H(z) has a frequency 
which is I times higher than the frequency X(z). Figure 5B shows the filter of Figure 5A 

decomposed into its polyphase components Hi : n(z), Hi : i_2(z), Hi : o(z). The polyphase 

components each receive the input signal X(z), and may therefore operate at the same 
frequency as the signal X(z) instead of a frequency I times higher. 

Unfortunately, even a symmetric filter may have asymmetric polyphase 
components, as is illustrated by the following example: 

Consider the filter h[n] having symmetric impulse response <a,b,c,d,d,c,b,a>. 
In case of multirate factor R=2, both its polyphase components h 2: o[n] and h 2: i[n] are non- 
symmetric: 

h2:o[n] = <a,c,d,b>, and 
h2:i[n] = <b,d,c,a>. 

It has been found by the inventors that the polyphase components of the 
multirate filter can be combined in (a,\\f) symmetric filter modules. The general structure of 
such a multirate filter according to the invention is illustrated in Figure 6A. 

The multirate filter according to the invention shown therein has: 

an input unit 1 0 for receiving an input signal Sin and for providing a plurality 
of intermediate signals IS1, IS2 in response to said input signal, 

a filter unit 20 coupled to the input unit 10, 

an output unit 30 coupled to the filter unit 20, for generating an output signal 

Sout. 

The filter unit 20 comprises at least a first and a second filter module 21, 22, 
having a transfer function Ho(z) and a transfer function Hi(z) respectively, which are 
mutually related according to the relations 
H 0 (z) = c 0 (H B (z) + M atV ,H B (z)) and 

H x 00 = c x (H B (z) -M ay ,H B (z)) . 

Therein M aJ9f H B (z) = az^H B \z' l ) 9 and H B \z) = £ h b \m]z' m . 
H B (z) is the z-transform of h b [m] . 
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By combining a non-symmetric filter Hb(z) with its mirrored version Ma iV 

H B (z), either a symmetric filter (in case of combining by addition) or an anti-symmetric filter 

(in case of combining by subtraction) is obtained. This is illustrated by the following 

example. When adding the two non-symmetric polyphase components 

h 2 :o[n] = <a,c,d,b>, and 
h 2: i[n] = <b,d,c,a>, 

the following symmetric transfer function is obtained: 
h[n] = <a+b,c+d,d+c,b+a> 

Subtracting these polyphase components from each other gives the following anti-symmetric 

transfer function. 

h'[n] = <a-b,c~d,d-c,b-a> 

The multirate filter comprises a combination unit 1 1 coupled to the filter 
modules 21, 22 for generating a first combination signal Ssum and a second combination 
signal Sdiff. The combination unit 1 1 may be a part of the input unit 10, or a part of the 
output unit 30. Alternatively both the input unit 10 and the output unit 30 may comprise one 
or more combination units. 

In an embodiment of the invention the basic transfer function Hb(z) is a 

polyphase component i.e. H B (z) = H R:r (z) = ^h[Rn + r]z" w 

n 

Figure 6B shows a first subembodiment of said embodiment. Therein the 
combination unit 1 1 is comprised in the input unit 10. The input unit 10 generates a first and 
a second intermediate signal IS1, IS2 from the input signal Sin. The first and the second 
intermediate signal IS1, IS2 are delayed by delay functions and z~* q respectively, wherein 
rq and sq are integer numbers selected from Z. After the delay the signals are decimated with 
a factor D, so that the intermediate signals have a sampling rate which is D lower than that of 
the input signal Sin. The combination unit 1 1 subsequently generates the first and the second 
combination signal Ssum, Sdiff from the pair of intermediate signals IS1, IS2. In the 
embodiment shown in Figure 6B the first combination signal Ssum is the sum of the 
intermediate signals IS1 and IS2 and the second combination signal Sdiff is the difference of 
the intermediate signals IS1 and IS2. The first filter module 21 filters the first combination 
signal Ssum and the second filter module 22 filters the second combination signal Sdiff. The 
output unit 30 subsequently combines the output signals of the filter modules 21 and 22 to 
generate the overall output signal Sout. The output unit 30 may additionally perform an 
interpolation with a factor I and introduce an extra delay z"* to the output signal Sout. 
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As is clear from the above, the transfer functions of the filter modules 21 and 
22 are symmetric. The transfer function Ho(z) of the first filter module is the sum of the basic 
transfer function H B (z) and its mirrored counterpart, and the transfer function Hi(z) of the 
second filter module is the difference of the basic transfer function H B (z) and its mirrored 
counterpart. Hence, both transfer functions Ho(z) and Hi(z) allow an efficient 
implementation, having a relatively small number of multipliers. On the one hand the filter 
modules can operate relatively slowly as they process the signal Ssum and Sdiff, which are 
derived from the decimated signals IS 1 and IS2. The optional interpolation stage only takes 
place after the filter unit 20. 

Figure 6C shows a second sub embodiment wherein the transfer function H B (z) 
is a polyphase component. In this case the first and the second filter module 121, 122 each 
filter a respective one of the intermediate signals IS1, IS2. The combination unit 131 is 
comprised in the output unit 130 and generates the first combination signal Ssum and the 
second combination signal Sdiff from output signals OS1, OS2 of the first and the second 
filter module. Optionally the combination signals are interpolated by interpolation factor I 
and delayed by a delay z" 1 " and z"*" respectively before they are recombined into the output 
signal. Likewise in this case the filter modules 121 and 122 are symmetric, and can be 
implemented with relatively small number of multiplier. Because the filter is operative after 
decimation and before interpolation, the filter can operate at a relatively low clock speed. 

The transfer functions Ho(z) and Hi(z) of the first and the second filter 
module are the sum of a basic transfer function H B (z) and its mirrored version Ma, v H B (z). In 
the embodiments shown in Figure 6B and 6C the basic transfer function is a polyphase 
component of the transfer function of the multirate filter. 

Figure 6D shows an embodiment wherein the basic transfer function H B (z) 
from which the transfer functions Ho(z) and Hi(z) are composed is derived from two 
polyphase components according to the relation: 
H B (z) = H R .^z) + z b H K . Ji (z). 

The combination unit 21 1 is comprised in the input unit 210, and generates the 
first and the second combination signal Ssum, Sdiff from a pair of the said intermediate 
signals IS1, IS2. The first filter module 221 filters the first combination signal Ssum and the 
second filter module 222 filters the second combination signal Sdiff. The filter unit 220 
further comprises a third filter module 223 with a transfer function H 2 (z) and a fourth filter 
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module 224 with a transfer function H 3 (z) which are mutually related according to the 
relations 

H 2 (z) = c 2 (H' B (z)-MH' B (z)) 
H 3 (z) = c,(H' B (z) + MH' B (z)). 
5 Also the further basic transfer function Hb'( z ) is derived from a combination 

of two polyphase components: 
H' B (z) = H R ^z)~z b H RA {z) 

The third 223 and fourth filter module 224 each filter a respective one of the 
first combination signal Ssum and the second combination signal Sdiff. The filter comprises a 
1 0 first further combination unit 23 1 for generating a first auxiliary signal AS 1 representing the 
sum of the output signals of the first filter module 221 and the second filter module 222, and 
a second auxiliary signal AS2 representing the sum of the output signals of the third filter 
module 222 and the fourth filter module 223. The filter further comprises a second further 
combination unit 232 for generating a further sum signal Ssum" and a further difference 

15 signal Sdiff" from the first and the second auxiliary signal AS1, AS2. These signals are 

i 

recombined into an output signal Sout after an optional interpolation with a factor I. The 

signals Ssum" and Sdiff' may further be delayed with an additional delay z~~ r ° p and z" ,oP . 

i 

As in the embodiments shown in Figures 6B and 6C the filter unit 220 of the 
multirate filter can on the one hand operate at a relatively low speed, as compared to a filter 

20 without polyphase decomposition. And on the other hand the modules 221 , 222, 223, 224 of 
the filter unit 220 have a symmetric transfer function, so that they can be implemented with a 
relatively small number of multipliers. 

Figure 7 A and 7B illustrate a further embodiment of the invention. Figure 7 A 
shows a multirate filter having a rational decimation factor 3/2, 

25 Figure 7B shows the inventive implementation of this filter. The filter unit 320 

comprises two groups of filtermodules. One group comprises the modules 232, 324, 325 and 
326, having transfer functions Eo(z), Ei(z), E2(z) and E 3 (z) respectively. The other group 
comprises the modules 321 and 322 having transfer functions Ba(z) and Es(z) respectively. 
The transfer functions for the first group are derived from the polyphase components of the 

30 original filter H(z) as follows. 

E 0 (z) = Uh 6 „{z) + H^ 2 (z) + H ^(z) + H 6i[ (z)) 
4 
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E x (z) = 1 {H M (z) + H M (z) - H a (z) - if frJ (z)) 

£ 2 (z) = i (tf M (z) - H & _ 2 (z) + tf 6 , (z) - H 6:1 (z)) 
4 

£ 3 (z) = 7(H«(*)-i5T«(*) - # «00 + &1 (z» . 
4 

The transfer functions for the second group are 
E,(z) = ^(H 6a (z) + H 6 ^z)) 

E 5 {z) = ±{H^(z)-H^{z)). 

It is noted that the skilled person will consider many alternatives within the 
scope of the invention as defined by the claims. For example the delay functions z 2 , z° and z 1 
in the input unit may be replaced by delay functions z 2 *, z' k and z l " k in combination with a 
delay function z k in the signal line providing the input signal X(z). It is noted that merely the 
relative delays generated between the signals provided by the delay functions in the input unit 
310 is of importance. 

It is also clear to the skilled person that several functions may be interchanged. 
For example, the decimate units indicated by the downpointing arrows with label 3 may be 
located at the outputs of the combination unit 311. Furthermore the delay units indicated with 
z° and z l may be included in the combination unit 31 1. In that case a first combination 
element of the combination unit 311 (the upper one in the drawing) generates the signal X(z) 
+ z ! X(z) = X(z) (1+z). The other combination element of said combination unit generates 
X(z)(l-z). 

The skilled person will furter realized that other operations can be 
interchanged for example the sequence of operations: 

interpolating with a factor I, 

delaying with z**, 

decimating with a factor D 
may be replaced by the sequence 

delaying with z aq 9 

decimating with a factor D 

interpolating with a factor I, 
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delaying with z**. 
Provided that D and I are coprime. 

Other examples are that multiplication operations can be interchanged with delay operations 
and with decimation or interpolation operations. 

A multirate filter according to the invention can be designed by carrying out 
the following steps: 

a. decomposing a symmetrical filter H into its polyphase components Hrm , 
...,Hr:r-i. 

This kind of decomposition is well known, and is described for example in the thesis 
"Efficiency in multirate and complex digital signal processing" by A.W.M. van den Enden, 
ISBN 90 6674 650 5. In particular reference is made to chapter 2.8 and to Appendix C 
thereof. 

The so obtained set of polyphase components may be replaced by a set of symmetrical 
modules as follows. Naturally those polyphase components which are already symmetrical 
need not be replaced. The following situations may occur. 

b. The set of polyphase components may comprise pairs of components Hr :i .i and 
H R:r2 which are asymmetrical, which share a common input, and which are related to each 
other according to the relation H R:r2 = M a<r (H R:rl ) . The parameters a,v|/ can be different for 

each pair. Such pairs of components Hr.ti and Hr^ can be replaced by 
a first module with transfer function H Q = H R:rl + H R:r2 , 
a second module with transfer function H x = H R:rl - H Rxr2 , wherein the first 

and the second module sharing the common input, and by 

a combination unit for generating a first and a second combination signal from 
the output signals of the first and the second module. 

c# The set of polyphase components may comprise pairs of components Hr*i and 

HR: r 2 which are asymmetrical, which share a common output, and which are related to each 
other according to the relation H R:r2 = M a ^ (H^ ) , wherein the parameters <x,i|/ can be 

different for each pair. 

Each such a pair of components Hr. t i and Hr. T 2 can be replaced by 

a first module with transfer function H 0 = H Rxl + H R:r2 , 

a second module with transfer function H % = H R:rl - H^ 2 , which modules 

share the common output, and by 

a combination unit for generating a first combination signal from a first and a 
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second intermediate signal, and providing said combination signal to the third 
module, and for generating a second combination signal from a first and a second 
intermediate signal, and providing said combination signal to the fourth module, 
d. Finally components may occur in quadruplets which are asymmetrical and 

which comprise 

a first pair Hr*i, Hr^ which are related by H R:r2 = M a t/f {H RtjrX ) 
a second pair Hr. t3 , Hr^ which are related by H M = M a%lff (H R:r2 ) 

wherein, 

the components Hr^i and Hr. t3 share a first common input, 
the components Hr^ and Hr. t4 share a second common input, 
the components Hr^ and Hr^ share a first common output, and 
the components Hr. t3 and Hrh4 share a second common output. 

Each such a quadruplet can be replaced by a quadruplet of symmetrical filter modules and 

three combination units. The quadruplet of filter modules comprises 

a first module having transfer function H 0 = H R ., 0 + H R:rl + H R:r2 + H Rxr3 , 
a second module having transfer function H x = H R:r0 - H Rlrl - H R:r2 + H R:rZ , 
a third module having transfer function H 2 = H R:r0 + H RirX - H Rt2 - H R:r3 , 
a fourth module having transfer function H 3 = H R:r0 - H R:rX + H R:r2 - H R:r3 . 
A first combination unit generates a first combination signal from input signals 

received at the first and the second common input, and provides said first combination signal 

to the first and the second module. It also generates a second combination signal from those 

input signals, and provides that to the third and the fourth module. 

A second combination unit generates a first auxiliary signal from output 

signals generated by the first and the third unit. It generates a second auxiliary signal from 

output signals generated by the second and the fourth unit 

A third combination unit generates a first and a second output signal from the 

first and the second auxiliary signal. 

It is noted that the situations b,c and d may occur in combinations. 



